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Introduction 

The  introduction  of  liquid  and  its  subsequent  at¬ 
omization  are  important  in  a  wide  variety  of  process 
and  subprocesses.  Understanding  atomization  is  there¬ 
fore  important  for  predicting  the  performance  of  many 
systems.  Atomization  is  generally  viewed  as  a  device¬ 
dependent  operation,  i.e.  atomization  quality  and  char¬ 
acter  is  considered  to  be  dependent  on  the  specific  type 
of  atomizer  utilized.  A  recent  review  of  the  mecha¬ 
nisms  effecting  film  atomization  [1]  brings  to  light  a 
device-independent  view,  however.  This  more  general 
viewpoint  has  been  applied  to  develop  an  atomization 
model  based  on  general  atomization  mechanisms. 
While  this  theory  does  not  apply  to  all  atomizers,  it  may 
be  utilized  to  predict  the  performance  of  many  different 
systems  provided  they  fall  within  certain  classes.  Spe¬ 
cifically,  the  theory  developed  here-in  applies  when 
droplets  are  produced  from  a  wall-bounded  film  and 
gas-phase  effects  are  important  due  to  high-speed  gas 
flows,  high  pressures  or  both.  As  will  be  shown,  the 
importance  of  the  gas  phase  can  be  estimated  via  a 
momentum-flux  ratio. 

The  specific  injector  motivating  this  study  is  a  gas- 
centered  swirl-coaxial  (GCSC)  injector.  A  schematic  of 
the  injector  is  shown  in  Fig.  1.  Liquid  is  injected  tan¬ 
gentially  along  the  outer  wall  of  the  atomizer.  This 
tangential  injection  causes  a  swirling  film  to  form  along 
the  wall.  High-speed,  nonswirling  gas  is  introduced 
axially  through  the  center  of  the  injector.  This  injector 
is  an  effective  atomizer  because  the  gas  flows  at  a  much 
higher  speed  than  the  liquid,  and,  therefore,  the  gas 
momentum  flux  is  larger  than  the  liquid  momentum 
flux.  These  injectors  share  characteristics  with  other 
atomizers  such  as  pressure-swirl  and  coaxial  air-blast 
atomizers,  but  unlike  these  other  injectors  they  do  not 
produce  a  conical  sheet  at  typical  operating  conditions. 
Instead,  atomization  occurs  from  a  wall-bounded  film 
inside  the  injector  cup  [2,  3],  More  information  on  this 
injector  as  well  as  some  experimental  and  numerical 
simulation  results  may  be  found  in  an  earlier  ILASS 
paper  [2],  Following  the  general  theoretical  develop¬ 
ment  an  example  is  given  using  this  atomizer. 

As  mentioned  above,  this  work  is  predicated  on  an 
earlier  examination  of  atomization  mechanisms  in  films 
where  several  likely  atomization  mechanisms  were  de¬ 
termined  [1],  Numerical  and  experimental  testing  of  a 
GCSC  injector  has  further  reinforced  the  identifications 
made  in  this  earlier  work  [2].  In  the  current  model 
droplets  are  produced  as  a  result  of  the  stripping  of 
various  disturbances.  These  disturbances  may  arise 
from  different  initiating  mechanisms  such  as  liquid  tur¬ 
bulence  or  hydrodynamic  instabilities.  This  paper  de¬ 
scribes  the  stripping  process  and  develops  a  theory  for 
calculating  the  rate  of  atomization,  atomization  effi¬ 
ciency  and  primary  droplet  size.  As  in  prior  simplified 
models  describing  film  breakup  in  cooling  tubes,  this 


theory  assumes  that  liquid  is  stripped  from  a  distur¬ 
bance  when  the  forces  on  a  section  of  the  disturbance 
are  balanced  [4,  5]. 

The  theoretical  development  is  presented  in  two 
subsections.  The  specifics  of  stripping  and  how  to 
quantify  it  are  given  first.  These  details  fonn  the  main 
backbone  of  the  theory  which  is  intended  to  be  as  gen¬ 
eral  as  possible:  it  applies  to  any  reasonably  sized  dis¬ 
turbance  in  a  flow  which  is,  on  the  average,  parallel  to 
the  bulk  liquid  interface.  The  second  subsection  pre¬ 
sents  the  procedure  for  calculating  atomization  effi¬ 
ciency.  Assumptions  and  procedures  for  calculating 
atomization  rates  and  primary  droplet  sizes  are  also 
given.  Following  the  presentation  of  the  general  theory, 
a  specific  example  is  given  for  a  GCSC  injector.  This 
example  serves  as  a  medium  for  discussing  various 
aspects  related  to  the  practical  application  of  the  theory. 


Figure  1.  Partial  schematic  of  the  experimental 
apparatus.  The  liquid  enters  through  the  tubes 
labeled  A;  gas  enters  at  B.  The  gas  inlet  extends 
beyond  what  is  pictured. 

Theoretical  Development 

Stripping 

In  the  simplest  sense  stripping  can  be  considered  to 
occur  as  one  of  three  different  modes — “pull”,  “push” 
and  “scoop”  (Fig.  2).  In  “pull”  stripping  the  lift  over  a 
curved  surface  causes  it  to  further  distort  and  eventually 
lifts  a  section  of  liquid  from  the  film  [6],  “Push”  strip¬ 
ping  is  the  result  of  drag  forces  where  a  section  of  a 
disturbance  is  dragged  from  the  film  [5].  “Scooping”  is 
the  result  of  gas-phase  structures,  such  as  a  recircula¬ 
tion  zone,  where  the  gas  motion  causes  a  distortion  in 
the  liquid — distorting  it  to  the  point  that  part  of  the  liq¬ 
uid  is  separated  from  the  bulk  [2].  “Scooping”  can  be 
thought  of  as  the  result  of  friction  instead  of  lift  or  drag. 


Figure  2.  The  three  modes  of  stripping  are 
illustrated  here:  (a)  “push”,  (b)  “pull”,  (c)  “scoop”. 


In  realistic  situations,  however,  stripping  is  an 
amalgamation  of  the  three  modes  since  all  liquid  protu¬ 
berances  are  subject  to  lift,  drag  and  other  forces.  De¬ 
spite  these  forces  acting  on  a  disturbance,  it  might  not 
undergo  stripping.  The  criterion  determining  if  strip¬ 
ping  occurs  is  based  on  a  simple  assumption:  if  the 
forces  over  a  section  of  the  disturbance  are  balanced 
then  that  section  will  separate  from  the  bulk  film.  This 
criterion  has  been  successfully  used  by  several  other 
researchers  to  calculate  atomization  rates  (due  to 
“push”  stripping  only)  in  cooling  tubes  with  multiphase 
annular  coflows  [4,  5]. 

Figure  3  shows  a  set  of  forces  on  a  generic  distur¬ 
bance.  The  general  direction  of  the  forces  is  indicated 
where  the  angles  are  only  approximate.  As  shown  in 
the  figure  lift,  drag  and  friction  forces  aid  atomization 
while  surface  tension,  centripetal  and  viscous  forces 
impede  the  separation.  Separation  occurs  at  some  angle 
to  the  radial  direction,  0S.  The  balance  of  forces  should, 
therefore,  be  considered  not  in  the  radial  or  axial  direc¬ 
tions  but  along  the  angle  6S. 

Lift  +  Drag( tan  9S )  +  Friction(cos  9 f  +  sin  6  f  tan  6S ) 

—  Centripetal  -  Surface  Tension^ cos  6a  +  sin  9a  tan  9S ) 

-  Viscous( cos  0 u  +  sin  0  tan  9S  )=  0 

(1) 

Since  the  trigonometric  terms  are  well-behaved  over  a 
wide  range  of  angles  (until  the  angles  approach  90°), 
the  exact  angle  need  not  be  specified:  each  trigonomet¬ 
ric  term  may  be  considered  a  constant  of  0(0.1)  to 
0(1).  Note  that  while  this  formulation  is  rather  general 
including  forces  which  may  not  occur  in  all  situations 
(e.g.,  friction  or  centripetal  forces),  it  is  slightly  simpli¬ 
fied  in  that  gravity  has  been  neglected.  The  gravity 
term  may  be  easily  incorporated,  but  care  must  be  taken 
that  the  term  is  included  in  the  correct  direction  as  it 
may  either  contribute  to  or  hinder  atomization. 

T  Lift 


Drag 


Figure  3:  A  depiction  of  the  forces  acting  on  the 
separating  portion  of  a  disturbance 

The  expressions  for  the  forces  need  not  be  exact; 
because  trigonometric  constants  are  already  present, 
adding  proportionality  constants  for  each  force  does  not 
appreciably  increase  the  formulation’s  complexity.  As 


a  result,  lift  and  drag  may  be  combined  into  a  single 
term  since  both  forces  are  proportional  to 

PgK^rel  (2) 

The  friction  term  is  similar,  but  the  specific  area  and 
velocity  involved  differ 

P  g  A  surface^  rel, retire 

Note  that  the  surface  area  is  not  that  of  the  entire  film  or 
disturbance,  but  the  area  in  contact  with  the  recircula¬ 
tion  zone  or  other  flow  feature  causing  the  friction 
force.  Similarly,  the  surface  area  in  the  viscous  force 
term  is  for  the  area  in  contact  with  the  distorting  flow 
field. 

PlASurfaceVrellS  (4) 

The  surface-tension  force  is  proportional  to 

/R  curv  (5) 

The  area  perpendicular  to  atomization  is  an  approxima¬ 

tion  and  is  rather  difficult  to  determine  since  the  atomi¬ 
zation  angle  is  not  resolved.  Consequently,  this  area 
will  be  simplified  to  either  the  annular  area  occupied  by 
the  disturbance  or  area  of  the  disturbance  seen  by  look¬ 
ing  along  the  axis  of  the  gas  flow.  Finally,  the  centripe¬ 
tal  acceleration  term  is  proportional  to 

PtVent®  centripetal  (6) 

with  the  centripetal  acceleration  given  by 

—  v2 

O' centripetal  *  local, 0 


/ rrot  .  If  included,  the  gravity  term  has 


the  same  form  as  the  centripetal  acceleration  term  with 
&  centripetal  replaced  by  the  gravitational  acceleration.  The 
various  areas,  entrained  volume  and  geometric  parame¬ 
ters  depend  on  the  disturbance  shape  and  atomizer  ge¬ 
ometry  (e.g.,  annular  or  flat). 

Utilizing  the  above  expressions  the  force-balance 
equation  contains  fluid  properties,  gas  and  liquid  ve¬ 
locities,  and  geometric  parameters  of  the  disturbance. 
Ideally,  all  of  these  are  known  except  one  geometric 
parameter — the  atomized  section  height,  5 — for  which 
the  equation  is  solved.  Practically,  however,  the  exact 
velocities  and  disturbance  shapes  are  rarely  known. 
Simulations  or  experiments  coidd  be  used  to  develop 
these  values  over  a  range  of  conditions,  but  this  ap¬ 
proach  is  unlikely  to  be  exhaustive  and  would  be  time 
and  cost  intensive.  It  is  recommended  instead  that  as¬ 
sumptions  be  made  in  order  to  approximate  these  not- 
fully-known  values.  In  particular,  the  shape  of  each 
disturbance  can  be  presumed  in  order  to  calculate  the 
various  geometric  parameters  found  in  the  force  terms; 
for  example,  a  disturbance  caused  by  hydrodynamic 
instabilities  could  be  assumed  to  have  a  sinusoidal  cross 
section.  The  relative  velocity  can  also  be  approximated 
as  the  gas  velocity  if  this  is  much  greater  than  the  liquid 
velocity.  Introducing  this  approximation  for  the  rela¬ 
tive  velocity  and  nondimensionalizing  with  p^inTgap , 
the  force-balance  equation  can  be  written  as 


®[Cz/D  Acs  j  Tgap  +  C f  (v Slrecirc  /Vg  )  surface  /  T  gap  )j 

“  RenVlocal,0  /{rrotT  gapV  ill  )_  C  aWe  /  A±atom  /{Rcun-T  gap 

-  CP  AsurfaceVg/{ReeVi„TgapS)=  0 

(7) 

where  <D  is  the  momentum-flux  ratio  of  the  gas  and 
liquid  at  inlet  conditions,  Wef  =  <j/(ptTgapvfn)  is  the 
liquid  Weber  number  and  the  Reynolds  number  of  the 
liquid  is  Re,,  =  P  P  ,,apv  in  j  p ,  .  This  equation  highlights 

the  importance  of  the  momentum-flux  ratio,  as  all  terms 
aiding  stripping  are  multiplied  by  it.  Consequently,  the 
assumption  that  the  gas  velocity  is  much  greater  than 
the  liquid  velocity  seems  justified  until  the  liquid  and 
gas  density  approach  one  another. 

General  Procedure  to  Calculate  Atomization  Rate 

The  first  step  in  determining  the  atomization  effi¬ 
ciency  and/or  average  droplet  diameter  is  to  develop  the 
necessary  experimental  measurements,  simulations  or 
assumptions  to  determine  what  type  or  types  of  distur¬ 
bances  are  present  on  the  film.  The  type  of  disturbance 
depends  on  the  operating  conditions  and  geometry  of 
the  atomizer  and  is  related  to  creation  mechanisms  as 
outlined  in  an  earlier  work  [1].  For  example,  liquid 
turbulence  creates  one  type  of  disturbance  (ligaments) 
while  hydrodynamic  instabilities  create  another 
(waves).  In  some  instances  a  single  type  of  disturbance 
will  be  present;  in  other  cases  numerous  types  will  exist 
on  the  same  film.  A  set  of  sample  questions  for  deter¬ 
mining  the  type  of  disturbances  present  is  presented  as 
Table  1. 

Once  disturbance  types  are  determined,  the  veloci¬ 
ties  and  geometric  parameters  appearing  in  the  force- 
balance  equation  must  be  established  for  each  distur¬ 
bance  type  present.  Again,  these  values  may  be  deter¬ 
mined  through  measurements,  simulations,  assumptions 
or  a  combination  thereof.  Generally,  common  sense 
and  prior  experiments  and  simulations  with  similar  at¬ 
omizers  or  conditions  will  allow  good  choices  of  dis¬ 
turbance  shapes.  Once  a  shape  is  chosen,  the  geometric 
parameter  can  be  determined  through  application  of 
geometry  and  simple  calculus.  The  determination  of 
the  radii  of  curvature  and  rotation  need  a  few  additional 
notes,  however.  Two  main  choices  present  themselves 
for  approximating  the  radius  of  rotation:  set  this  radius 
equal  to  the  average  radius  over  the  atomized  section  of 
the  disturbance  (simple)  or  set  it  equal  to  the  radial 
value  of  the  cross-section’s  centroid  (accurate).  The 
radius  of  curvature  is  more  complex  and  an  obvious 
assumption  does  not  present  itself.  For  ease  and  be¬ 
cause  of  its  successful  use  in  earlier  force-balance  for¬ 
mulations  [5],  the  value  at  the  crest  of  the  disturbance  is 
recommended  as  an  approximation  for  the  value  at  at¬ 
omization. 


Can  aerodynamic  effects  clearly  be  neglected?  <D<Q(T)* 
YES — this  mechanism  does  not  apply 
NO — move  on  to  the  disturbance  type  questions 
Stripping  from  liquid  turbulence  disturbances  possible? 
(must  answer  yes  to  all  below) 

Is  the  liquid  turbulent? 

Do  the  eddies  have  sufficient  energy  to  deform  the 
surface  (Wes;My>0(l))r! 

Is  the  outlet  length  sufficient  for  the  eddy  to  fonn? 
Stripping  from  hydrodynamic  instabilities  waves  possi¬ 
ble?  (must  answer  yes  to  all) 

Is  the  flow  unstable  or  transiently  unstable? 

Is  the  injector  of  sufficient  length  for  a  wave  to  grow 
to  a  height  where  it  could  be  atomized? 

Are  disturbances  due  to  a  recirculation  zone  possible? 
(must  answer  yes  to  all) 

Is  there  a  potential  recirculation  zone  due  to  the 
geometry  of  the  injector  (step,  steep  slope,  other 
area  of  possible  separation)? 

Is  the  recirculating  flow  energetic  enough? 

(d>(vre/ve)2>0(l)*) _ 

Are  disturbances  due  to  gas-phase  turbulence  possible? 
(must  answer  yes  to  all) 

Is  the  gas  turbulent? 

Is  the  flow  energetic  enough  to  be  important? 

(d>  (Vm,s/Ve)2>0(l)*) _ 

*in  some  cases  the  criterion  will  be  different 

Table  1 :  Sample  questions  for  determining  possible 
disturbance  types.  Note  that  existence  of  a  distur¬ 
bance  does  not  guarantee  stripping  of  the  distur¬ 
bance.  For  example,  in  strong  gas  flows  ligaments 
produces  by  liquid  turbulence  may  topple  before 
they  can  be  stripped. 

Velocity  determinations  tend  to  be  atomizer  specific 
and  require  either  a  large  body  of  experimen¬ 
tal/simulation  work  or  many  assumptions.  Flowever, 
progress  can  be  made  with  straightforward  approxima¬ 
tions.  For  example,  the  local  tangential  velocity  of  the 
liquid  may  be  considered  to  vary  linearly  with  distance. 
The  gas  velocity  may  be  equated  with  that  of  the  in¬ 
coming  gas,  related  to  the  distance  via  correlations  for 
flow  over  wavy  surfaces  [5]  or  estimated  from  single¬ 
phase  simulations  of  the  atomizer  [7].  Experimentally 
measured  values  could  also  be  used,  but,  unless  some 
general  correlation  can  be  determined  from  a  limited 
number  of  experiments,  this  approach  would  be  very 
time  consuming  and  costly. 

Once  the  velocity  and  geometric  values  are  estab¬ 
lished,  Eq.  7  may  be  solved  for  the  height  of  the  distur¬ 
bance  lost  to  atomization,  8.  This  lost  height  along  with 
the  shape  approximation  gives  a  value  for  the  entrained 
volume.  The  calculated  entrained  volume  is  that  lost  by 
an  individual  disturbance.  Several  protrusions  may  be 
present  on  the  fdm,  however,  so  the  number  of  each 
type  of  disturbance  must  also  be  ascertained.  This 


number  is  related  to  the  film  length,  volume  or  surface 
area  of  the  liquid  and  the  disturbance  type.  Axial  film 
velocities  and  atomizer  geometry  may  also  play  a  role. 

The  total  atomized  volume  is  the  sum  of  the  liquid 
atomized  from  all  of  the  disturbances.  The  simplest 
approach  is  to  estimate  the  stripping  as  a  global,  i.e. 
time-independent,  value.  In  most  cases,  however,  the 
stripping  is  cyclic  since  a  disturbance  forms,  grows, 
loses  mass,  grows  and  loses  mass  again.  An  unsteady 
formulation,  while  more  complex,  allows  this  shedding 
frequency  to  be  estimated.  A  time-dependent  approach 
often  makes  the  calculation  of  the  number  of  distur¬ 
bances  on  the  film  more  straight-forward,  particularly 
in  the  situation  where  disturbances  grow  as  they  travel 
downstream.  Waves  created  by  hydrodynamic  insta¬ 
bilities  are  one  example  of  such  a  disturbance.  The 
growth-rate  of  a  disturbance  varies  depending  on  the 
cause  of  the  disturbance.  The  time  necessary  for  a  sec¬ 
tion  of  liquid  to  be  stripped  from  the  disturbance  (once 
the  disturbance  reaches  a  particular  height)  is  assumed 
to  be  the  time  for  the  gas  to  flow  over  the  atomizing 
portion  of  the  disturbance.  Using  the  knowledge  of  the 
number  of  disturbances  on  the  film,  their  growth  rate 
and  an  atomization  time  allows  the  calculation  of  a 
volumetric  atomization  rate — 

atom  ^ disturb^ en  / atom  (^) 

disturb  _  type 

The  atomization  efficiency  is  found  by  comparing 
the  volumetric  flow  rate  of  the  incoming  liquid  with  the 
atomization  rate  from  all  disturbances.  In  most  cases,  a 
film  length  is  needed  to  arrive  at  the  total  number  of 
disturbances.  To  get  the  correct  atomization  rate,  the 
correct  film  length  must  be  used.  This  film  length  is 
unknown  at  the  start  of  the  calculations  and  must  be 
estimated.  If  the  atomization  rate  calculated  from  the 
estimated  film  length  is  not  equal  to  the  inlet  flow  rate 
then  the  film  length  is  incorrect  and  a  revised  estimate 
for  the  film  length  should  be  used  to  get  a  new  atomiza¬ 
tion  rate.  This  iteration  process  continues  until  the  two 
flow  rates  are  equal.  Note,  however,  that  the  film 
length  cannot  exceed  the  atomizer  length  as  this  theory 
is  built  on  the  assertion  that  atomization  occurs  from  a 
film  only.  Once  the  correct  film  length  has  been  deter¬ 
mined  then  the  atomization  efficiency  can  be  calculated 
as  [7] 

atom  —  tnj  ~  ^  film  )/^Jinj  00 

Atomization  efficiency  is  only  one  part  of  atomizer 
performance;  the  droplet  size  is  another  important  pa¬ 
rameter.  The  force  balance  results  in  the  calculation  of 
an  atomized  volume  from  each  disturbance.  Stripping 
is  assumed  to  occur  across  the  entire  disturbance  which 
may  be  the  length  of  the  film  or  some  shorter  distance. 
Many  types  of  disturbances  will,  therefore,  produce  a 
ligament  when  stripping  occurs;  the  length  of  these 
ligaments  equals  the  disturbance  length.  For  example, 


an  annular  disturbance  would  produce  a  toroid-shaped 
ligament.  Using  the  ligament  length  and  the  calculated 
entrained  volume,  the  diameter  of  the  ligament  is 


(Note  that  for  if  a  toms  is  produced  LUg  is  the  circum¬ 
ference  at  the  center-point  of  the  atomized  section.) 
Assume  that  this  ligament  breaks  up  via  the  Rayleigh 
mechanism,  so  that  it  yields  droplets  with  a  diameter 
-1.89  times  the  diameter  of  the  ligament  [8].  The  num¬ 
ber  of  droplets  produced  by  this  ligament  is  roughly 

N drop  =  /[^(l  -89rf/ig  j3  J  (11) 

This  droplet-number  equation  may  be  used  in  all  cases, 
even  when  a  single  droplet  is  produced;  in  that  case  the 
number  of  droplets  calculated  will  be  one  or  less. 
When  the  number  of  droplets  is  greater  than  one  but  not 
an  integer,  the  ligament  should  be  assumed  to  break 
into  somewhat  differently  sized  droplets  in  order  to 
generate  an  integer  number.  In  other  words,  the  num¬ 
ber  of  droplets  should  be  rounded  to  the  nearest  integer. 
The  droplet  diameter  is  then 

Jf/3  (i2) 

where  Njrop  is  the  integer  value.  As  with  the  entrained 
volume  calculation,  the  entire  droplet  size  procedure  is 
applied  to  each  type  of  disturbance.  The  average  drop¬ 
let  diameter  produced  by  primary  atomization  is  a 
weighted  average  of  the  calculated  droplet  sizes 

d avg  ^ dropd drop  I  I ]Ndrop  (13) 

disturb  _  type  /  disturb  _  type 

where  the  summations  are  taken  over  all  disturbance 
types.  A  summary  of  the  entire  process  is  given  in  Ta¬ 
ble  2. 

Results  and  Discussion 

As  pointed  out  above,  the  application  of  this  theory 
requires  information  that  is  not  readily  known  a  priori 
and  may  be  specific  to  the  atomizer.  In  order  to  dem¬ 
onstrate  this  theory  in  action,  therefore,  an  example  is 
given  of  a  specific  atomizer — a  GCSC  injector.  The 
geometry  and  general  operation  of  a  GCSC  injector  are 
described  in  the  introduction  and  shown  in  Fig.  1.  Due 
to  space  constraints  and  its  illustrative  nature,  this  ex¬ 
ample  case  is  simplified  by  considering  only  a  small 
range  of  operating  conditions  where  a  single  type  of 
disturbance,  one  caused  by  gas-phase  turbulence,  is 
important.  In  this  regime  eddies  in  the  gas  are  so  ener¬ 
getic  that  they  distort  the  surface,  producing  protru¬ 
sions.  Example  operating  conditions  in  this  range  as 
well  as  some  nondimensional  parameters  are  given  in 
Table  3. 

The  first  information  that  depends  on  the  specific  at¬ 
omizer  (and  is  often  not  fully  known)  are  the  relative 
and  liquid-tangential  velocities.  The  gas  velocity, 
which  is  purely  axial  in  this  example,  is  much  greater 


Table  2:  A  summary  of  the  basic  process  for  determining  atomization  efficiency  and  average  droplet  diameter. 


than  the  liquid  velocity.  Therefore,  assume  the  bulk 
relative  velocity  is  also  purely  axial  and  equal  to  the  gas 
velocity.  Earlier  work  on  a  similar  atomizer  suggested 
that  the  appropriate  gas  velocity  is  not  its  bulk  velocity 
but  its  value  near  the  interface  which  can  be  approxi¬ 
mated  from  a  single-phase  numerical  simulation  of  the 
injector  [7].  The  relative  velocity  of  recirculating  gas, 
important  in  the  friction  term,  arises  from  turbulent 
eddies  and  can  be  estimated  from  turbulent  statistics 
such  as  the  root-mean-square  velocity  fluctuation  in  the 
axial  direction.  Multiphase,  VOF  simulations  (see  Fig. 
4)  suggest  that  the  liquid-tangential  velocity  falls 
sharply  following  liquid  injection  and  that  the  ratio  of 
tangential  to  inlet  velocity  is  0(1);  an  initial  approxima¬ 
tion  might  assume  Viocai:g=Cgvin.  From  simulations  a 
more  accurate  functional  dependence  between  the  two 
velocities  could  be  developed;  however,  these  simula¬ 
tions  are  quite  time  consuming  and  given  the  other  as¬ 
sumptions,  unlikely  to  appreciably  improve  accuracy. 

Assume  that  the  eddies  contacting  the  liquid  sur¬ 
face  are  spherical  or  cylindrical  in  shape.  Turbulence 
literature  [9]  suggests  the  shape  may  be  quite  a  bit  more 
complex;  yet,  atomization  literature  studying  atomiza¬ 
tion  due  to  liquid  turbulence  has  been  quite  successful 


Gap  Thickness  (in) 

0.065 

Outlet  Radius  (in) 

0.375 

Liquid  Mass  Flow  (lb/s) 

0.101 

Gas  Mass  Flow  (lb/s) 

0.135 

Liquid  inlet  velocity  (ft/s) 

10-30 

Gas  inlet  velocity  (ft/s) 

1000  + 

d) 

3-5 

We, 

0(1  O'3) 

_ 

O(104) 

Table  3:  Operating  conditions  for  the 
GCSC  injector.  Working  fluids  are 
gaseous  nitrogen  and  water;  test  con¬ 
ducted  at  STP. 

using  simplified  eddy  shape  and  behavior  [8].  The  area 
and  volume  calculations  are  most  straightforward  if  the 
eddy  is  assumed  to  be  cylindrical;  a  rough  length  esti¬ 
mate  may  be  obtained  from  tangential  velocity  fluctua¬ 
tions  or  by  assuming  the  eddy  is  somewhat  equivalent 
to  a  sphere  and,  therefore,  has  a  length  equal  to  half  its 
cross-sectional  circumference.  This  rotating  eddy  pulls 
liquid  up  along  its  outer  edge  in  a  manner  similar  to  a 
rotating  dip  coater.  The  resulting  disturbance  is  illus- 


trated  in  Fig.  2c;  as  seen  in  the  figure,  determining  the 
relevant  areas  involves  several  unknown  parameters 
such  as  Reddy,  a  distance  from  the  center  of  the  eddy  to 
the  nominal  film  surface  (y)  and  a  radius  of  the  outer 
portion  of  the  disturbance  ( Rout ).  The  eddy  radius  can 
be  determined  from  turbulence  statistics  and  the  other 
parameters  can  be  calculated  by  applying  simple,  logi¬ 
cal  conditions.  For  example,  the  volume  of  fluid  dis¬ 
placed  above  the  original  surface  is  equal  to  the  volume 
of  the  divot  (the  void  below  the  undisturbed  film  sur¬ 
face).  Secondly,  assume  that  the  liquid  is  drawn  to  the 
top  of  the  eddy  at  the  time  of  atomization  then 
h~Reddy+y ■  (This  conditions  means  the  center  of  the 
outer  circle  must  lie  on  the  same  vertical  line  as  the 
circle  describing  the  eddy.)  To  fully  describe  the  ge¬ 
ometry  one  final  condition  is  necessary  or  the  force- 
balance  equation  may  be  solved  to  get  Jas  a  function  of 
one  of  the  geometric  parameters,  such  as  h  or  Rout.  Cal¬ 
culations  using  simplified  geometries,  the  values  cited 
in  Table  3,  vnns  equal  to  30%  of  the  bulk  gas  velocity, 
several  values  of  Reddy  starting  from  the  minimum  dis¬ 
cussed  in  the  following  paragraphs  and  a  range  of  h 
indicate  that  h  must  be  relatively  large,  larger  than  Rouh 
before  atomization  occurs.  If  a  definitive  answer  is 
desired  rather  than  a  possible  range  as  h  varies,  a  condi¬ 
tion  equating  the  kinetic  energy  of  the  eddy  and  the 
surface  energy  of  the  overall  disturbed  shape  is  sug¬ 
gested. 
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Figure  4:  The  tangential  velocity  plot  from  an  axi- 
symmetric,  two-phase  simulation  of  the  injector. 

The  black  line  represents  the  approximate  location 
of  the  interface.  Note  that  this  figure  is  for  a  gas 
velocity  around  100  m/s  only. 

The  force-balance  equation  can  now  be  formulated 
in  terms  of  two  still-unknown  parameters — 8  and  Reddy- 
In  many  cases  the  problem  may  be  simplified  by  ne¬ 
glecting  certain  forces  that  are  much  smaller  than  the 
others.  Forces  can  be  calculated  with  8  values  near  0 
and  near  h  using  the  operating  conditions  and  either  the 
above  complex  geometry  or  a  simplified  but  similar  one 
(e.g.,  a  triangle  of  height  8  and  a  width  similar  to  the 
atomized  portion  of  the  disturbance).  If  any  force  is 
small  over  the  range  of  8  it  may  be  neglected;  other 
forces  may  be  neglected  if  small  at  either  limit  and  then 
rechecked  after  a  value  of  8  is  calculated.  If  all  of  the 
positive  terms  are  negligible  then  stripping  does  not 
occur.  Somewhat  surprisingly,  the  sample  calculations 


mentioned  above  (parameters  as  in  Table  3,  v,.ms=0.30vg 
and  a  range  of  h  and  Reddy)  indicate  that  the  surface 
forces  are  negligible  in  most  cases,  generally  being  two 
or  more  orders  of  magnitude  smaller  than  viscous  and 
centripetal  forces. 

The  amount  of  liquid  atomized  by  a  single  distur¬ 
bance  may  now  be  calculated  if  the  general  size  of  the 
eddy  is  known.  Since  this  calculation  is  for  a  single 
eddy,  the  number  of  disturbances  must  also  be  deter¬ 
mined  in  order  to  calculate  the  total  atomized  volume. 
Clearly,  the  size  and  number  of  eddies  depends  on  the 
flow  conditions  and  type  of  turbulence.  In  this  case 
turbulence  is  a  result  of  developed  pipe  flow.  In  order 
to  distort  the  interface  the  energy  of  the  turbulent  eddy 
must  exceed  the  surface  energy  of  the  film,  deddy 
~  olpveddy',  this  follows  the  theoretical  development 
predicting  disturbance  sizes  due  to  liquid-phase  turbu¬ 
lence  [8].  If  the  eddies  responsible  for  disturbance  for¬ 
mation  lie  in  the  inertial  range,  as  they  do  for  liquid 
turbulence,  then  vedd}~urms{ deddyJd /) 1 3  where  the  root- 
mean-square  velocity  is  in  the  cross-stream  direction 
and  d[  is  the  integral  length  scale,  dHl 8  [8].  Eddies  re¬ 
sponsible  for  atomization  therefore  meet  the  following 

Criterion  deddy/ dfi  Cedd^S^I pd} jVeddy^)  eddy  W^d / /, eddy  , 

at  a  minimum,  with  the  proportionality  constant  0(1). 
Starting  with  this  minimum  eddy  size,  with  a  constant 
of  unity,  results  in  a  primary  droplet  size  around  20-150 
pm.  Droplet  size  measurements  with  similar  atomizers 
[7]  found  SMDs  less  than  75  pm  to  as  fine  as  5  pm. 
Given  that  these  measurements  were  made  downstream 
of  the  atomizer  exit,  after  secondary  breakup  occurred, 
the  theory’s  findings  seem  reasonable.  The  percentage 
of  eddies  meeting  the  criteria  for  producing  distortions 
as  well  as  the  average  distortion-producing  eddy  size 
can  be  estimated  from  the  turbulence  spectrum  of  the 
atomizer  or  a  similar  configuration,  such  as  a  rough 
pipe. 

Of  the  eddies  capable  of  deforming  the  surface 
only  a  fraction  actually  impact  the  film;  this  fraction  is 
related  to  the  amount  and  profile  of  liquid  in  the  injec¬ 
tor.  Imagine  that  the  profile  of  the  film  is  roughly  lin¬ 
ear  with  a  maximum  at  the  lip  edge  where  the  film 
thickness  is  the  gap  thickness,  TiOCarTgap(l-z/Lf, i„,).  The 
film  volume  is 

Vfiim  ~  KL  fUmT gap  (3r0  —  rgap  )/3  (17) 

Assume  that  the  ratio  of  the  film  volume  to  total  vol¬ 
ume  is  equivalent  to  the  fraction  of  existing  eddies  to 
eddies  which  impact  the  film.  The  total  number  of  dis¬ 
turbances  formed  is 

N  -  1  Lfilm  T  gap&o  ~  T  gap)  (0/  ,r  j  /ion 

V  disturb  ,  2  \/0 capable 1  ’  eddy /  v10/ 

j  L .  .  r 

inj  o 

Unfortunately,  the  overall  film  length  is  not  initially 
known  and  must  be  guessed  in  order  to  calculate  these 
parameters.  Simulations,  experiments  or  prior  calcula- 


tions  are  suggestions  for  developing  an  initial  guess  of 
the  film  length.  In  the  absence  of  other  information  the 
injector  length  is  a  good  starting  point. 

As  for  atomization  times,  suppose  that  at  any  instant 
the  total  number  of  impacting  eddies,  Ndisturb,  are  in  con¬ 
tact  with  the  film.  These  eddies  require  a  finite  time  to 
distort  the  film  and  achieve  atomization;  the  next  batch 
of  eddies  then  contact  the  film  after  some  lag.  In  actu¬ 
ality,  this  type  of  disturbance  formation  and  atomization 
are  continual  processes  with  no  set  atomization  fre¬ 
quency;  however,  this  simplification  makes  the  calcula¬ 
tion  of  an  atomization  time  more  straightforward. 
Forming  the  distortion  requires  a  time  of 
dReddJvg, retire — the  time  for  the  circulating  flow  within 
the  eddy  to  rotate  up  along  the  disturbance  height.  The 
atomization  event  time  is  that  needed  for  the  circulating 
flow  to  move  across  the  atomizing  portion  of  the  dis¬ 
turbance.  The  lag  time  between  batches  of  eddies  is 
needed  to  account  for  the  reality  that  all  of  the  eddies 
do  not  contact  the  film  at  once.  It  may  be  approximated 
as  an  average  traveling  time  for  an  eddy  to  move  from 
its  position  in  the  gas  to  the  liquid  surface,  (r,J2)/urms. 
Overall,  then,  the  process  is  approximated  as  Ndisturb 
atomizations  occurring  every 

@R-eddy/v g, retire^" g, retire^- ^ d (^Urms)  Seconds.  Combining 
the  above  information 

Qatom  ^ disturb  F ?d t  (19) 

Where  OReddJVg, retire^ d/Vg,recirc~^~^ <J (2ll rms)' 

The  calculated  atomization  rate  is  compared  with  the 
volumetric  flow  rate  of  the  liquid  into  the  injector.  If 
these  values  are  not  equal  then  the  film  length  is  incor¬ 
rect  and  a  new  estimate  is  needed.  This  estimate  may 
be  based  on  the  atomization  rate  directly,  or  it  may  be 
calculated  by  considering  the  difference  in  volumetric 
flow  rates.  The  first  choice  is  slightly  less  complex,  but 
tends  to  overestimate  the  difference  between  the  lengths 
and  require  more  iterations.  To  calculate  a  film  length 
based  on  the  difference  in  volumetric  flow  rates,  let 
AQtJIoKr=AV.  The  flow  time  can  be  represented  by 
L film  I V local  Z  where  the  bar  represents  an  average 

value  over  the  film’s  length.  The  average  axial  velocity 
is  the  integral  of  the  local  axial  velocity  over  the  film 
length.  Setting  the  volumetric  flow  rate  through  the  r-9 
plane  as  a  constant  equal  to  the  inlet  flow  rate  and  fur¬ 
ther  assuming  the  atomization  is  evenly  distributed  over 
the  film  gives  an  expression  for  the  local  velocity: 

V local,:  =  Qin  / ^local  (^0) 

Afocai  is  just  the  annular  area  from  the  pipe  wall  to 
Tiocai-  The  volume  difference  is  Eq.  17  with  Lfltm  re¬ 
placed  by  ALj,/m.  The  difference  in  flow  rate,  AQ,  is  just 
Qin-Qatom •  Combining  all  of  these  pieces  allows  an 
equation  for  ALjum  to  be  developed. 

Once  the  calculated  atomization  and  inlet  flow  rates 
are  equal  then  the  film  length  has  been  correctly  calcu¬ 


lated.  The  atomization  efficiency  is  then  given  by  Eq. 
9.  Since  only  one  disturbance  type  is  considered  the 
average  droplet  diameter  due  to  primary  atomization 
may  be  found  by  solving  the  force-balance  equation  for 
d  and  using  the  entrained  volume  formula  along  with  a 
ligament  length  of  KcleddyJ2. 

Following  the  above  outline,  if  the  turbulence 
parameters  are  known  or  may  be  determined  (say  from 
prior  experiments/DNS  of  fully  developed  pipe  flow) 
then  the  atomization  efficiency  and  average  primary 
droplet  size  may  be  determined.  As  indicated  above, 
sample  calculations  suggest  primary  droplet  sized  of 
20-150  pm  for  the  atomizer  examined  here.  The  size 
and  amount  of  disturbance  atomized  suggests  that  up  to 
40%  of  the  liquid’s  surface  must  be  covered  with  atom¬ 
izing  disturbances  to  fully  atomize  the  film  over  the 
length  of  the  atomizer.  Since  the  ratio  of  film  to  total 
volume  is  only  25%  this  suggests  that  the  primary  drop¬ 
let  size  is  likely  to  be  in  the  upper  value  of  the  calcu¬ 
lated  range  and  that  eddies  larger  than  the  minimum 
necessary  size  likely  play  an  important  role  in  atomiza¬ 
tion. 

Conclusions 

A  theory  has  been  developed  which  uses  an  under¬ 
standing  of  general  atomization  mechanisms  to  calcu¬ 
late  atomization  properties  of  devices  utilizing  films 
and  having  strong  gas-phase  involvement.  A  general 
outline  of  the  process  was  given  as  well  as  an  example 
from  a  specific  atomizer  operating  in  a  restricted  re¬ 
gime.  The  theory  assumes  that  atomization  occurs 
through  a  process  of  surface  disturbance  creation  fol¬ 
lowed  by  the  gas-phase-initiated  stripping  of  a  portion 
of  these  disturbances.  The  model  is  applicable  to  any 
number  of  disturbance  creation  mechanisms;  only  one, 
that  due  to  gas-phase  turbulence,  was  considered  in 
detail  here.  Stripping  occurs  when  the  forces  acting  on 
a  small  section  of  a  disturbance  are  balanced.  This 
small  section  separates  from  the  bulk  fluid,  forms  a 
ligament  and,  eventually,  droplets.  Application  of  the 
theory  requires  several  assumptions,  specifically  deter¬ 
mining  velocities  and  disturbance  shapes.  Disturbance 
shapes  can  often  be  estimated  rather  accurately  while 
the  velocity  estimations  remain  a  large  source  of  uncer¬ 
tainty.  The  uncertainty  may  be  mitigated  through  ex¬ 
perimental  studies  or  by  the  use  of  a  limited  number  of 
simulations  to  determine  suitable  velocity  relations. 
Overall,  however,  the  theory  is  powerful  in  its  ability  to 
offer  predictions  for  a  wide  range  of  atomizers  operat¬ 
ing  over  a  wide  range  of  conditions.  Initial  calculations 
are  promising,  predicting  primary  droplet  diameters  in 
an  expected  range.  The  model  also  fills  a  void  in  that  it 
deals  with  the  film,  or  wall-bounded  liquid,  configura¬ 
tion  which  has  received  relatively  little  attention  and, 
therefore,  has  few  empirical  correlations  or  analytical 
prediction  methods. 


A  large  amount  of  future  work  is  planned.  Much 
work  is  needed  to  validate  this  model  and  find  suitable 
ranges  for  the  constants  contained  within  it.  Work  is 
also  planned  to  further  elucidate  the  effects  of  various 
assumptions  and  to  allow  recommendations  of  suitable 
assumptions  over  a  wide  range  of  conditions  and  ge¬ 
ometries.  The  theory’s  sensitivity  to  various  operating 
parameter  will  also  be  explored. 
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Nomenclature 

a  acceleration 

A  Area 

C  constant 

d  diameter 

h  overall  height  of  a  disturbance 

L  length 

N  number  of  things 

Q  volumetric  flow  rate 

r,R  radius 

Re  Reynolds  number 

t  time 

u  radial  velocity 

v  velocity 

V  volume 

y  distance,  see  Fig.  5 

z  axial  coordinate 

We  Web  er  numb  er 

5  height  of  the  atomized  portion  of  the  distur¬ 

bance 

r|  efficiency 

9  an  angle 

|i  viscosity 

p  density 

ct  surface  tension 

t  thickness 

<D  momentum-flux  ratio 

Subscripts 

_L  atom  perpendicular  to  atomization 

atom  atomization 

avg  average 

cs  cross-section 

curv  curvature 

disturb  disturbance 

drop  droplet 

eddy  turbulent  eddy 

en  entrained 

f  friction 

g  gas 

gap  of  the  gap  prior  to  where  the  liquid  and  gas 
meet 


in  inlet  value  of  liquid 

inj  injector 

1  liquid 

L/D  Lift  and  Drag 

lig  ligament 

recirc  recirculation 

rel  relative 

rot  rotation 

s  at  which  stripping  occurs 

o  outlet  of  the  injector 

z  in  the  axial  direction 

recirc  recirculation 

rel  relative 

p  viscous 

0  in  the  angular  direction 

ct  surface  tension 
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